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Abstract 

We analyse statistical properties of the normal forms of random braids 
generated by two different methods, namely generating random words of 
given length in the generators, respectively generating random braids with 
uniform distribution on the set of braids of given length. 

Except for an initial and a final region whose lengths are uniformly 
bounded, the distributions of the factors of the normal form of sufficiently 
long random braids depend neither on the position in the normal form nor 
on the lengths of the random braids. Moreover, when multiplying a braid 
on the right, the expected number of factors in its normal form that are 
modified is uniformly bounded; the latter property yields an algorithm for 
computing normal forms that has linear expected run time. 



1 Introduction 

Explicit computations play an increasingly important role in most areas of al- 
gebra; the study of braids is no exception. In many situations, computations 
with braids involve choosing braids at random: Some algorithms explicitly re- 
quire a random braid to be generated; this is the case, for instance, in crypto- 
graphic protocols based on the braid group |AAG99[ IKLC+OO] . At other times, 
a large collection of typical examples is to be generated; this is usually the case 
in computational experiments supporting theoretical research. 

As _B„, the group of braids on n strands, is infinite, choosing braids at random 
is not a trivial task. There are various natural ways of choosing elements of _B„ 
at random, and different approaches will yield different probability distribu- 
tions on Bn. For both, computational experiments and applications (especially 
applications in cryptography), it is important to understand the statistical prob- 
abilities of samples of random elements generated using a particular method. 

We consider in the following the braid monoid defined by the presentation 

CT,(7j=cr^CT, {l<i<j<n)\^ 
(1) = ( 0-1, . . . ,cr„_i l^ ^ ■ ^ 1^ / • 
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As the relations of -B+ are homogeneous, the number of generators occurring in 
any expression of x G is well-defined; we call this number the length \x\ of x. 
We can then fix a non-negative integer k and generate an element x G of 
length k. More specifically, there are two possibilities: 

(A) For i = 1, . . . , fc independently choose ai A — {cti, . . . , (Tn-i} with a 
uniform probability distribution on A, or equivalently, consider the set 
A*\k of all words of length fc over the alphabet A, and choose an element 
of A*\k at random with a uniform probability distribution on this set. 

(B) Consider the set B^\k — {x ^ B^ : \x\ — fc} and choose an element of 
B^\k at random with a uniform probability distribution on this set. 

We will refer to (A) as generating uniformly random words, and to (B) as gener- 
ating uniformly random braids. We will write Wordfc, respectively URB^, for the 
corresponding probability measures on B^ . Since the number of different words 
in A*\k that represent the same element x of B^\k depends on x, generating 
uniformly random words results in a distribution of braids which is very far from 
being uniform on B^\k. Generating uniformly random braids is not easy; an 
algorithm whose time-complexity respectively space-complexity is polynomial 
in both n and fc was given in |GGM13] . 

In this paper we compare the generation of uniformly random braids to the gen- 
eration of uniformly random words regarding some properties of the generated 
samples of braids. The Garside normal form defines a canonical way of express- 
ing a braid as a sequence of permutations, so a probability distribution on the 
braid group induces a sequence of probability distributions on the symmetric 
group. We are particularly interested in how the resulting distributions on the 
symmetric group depend on the position in this sequence. 

The structure of the paper is as follows: [section 2l recalls the Garside normal 
form; experts may skip this section, [section Sl contains our analysis of the normal 
forms of random braids. In lsubsection 3.ll we show that there is a "stabilisation" 
occurring in the normal forms of long random braids in the sense that for suf- 
ficiently long braids, the distributions on the symmetric group induced by the 
factors of the normal form depend neither on the position in the normal form 
nor on the lengths of the random braids, except for an initial and a final region 
whose lengths are uniformly bounded. In lsubsection 3.2l we give an explanation 
for this stabilisation phenomenon by demonstrating that the expected number 
of factors of the normal form of a braid that are modified when multiplying the 
braid on the right is uniformly bounded. Finally, in lsubsection 3.3l we extend 
our analysis to general Garside groups and establish a criterion for deciding 
whether the above mentioned phenomena occur in a given Garside group. 

2 Background 

This section contains a brief summary of the main notions referred to in the 
paper. Specifically, we will recall Garside monoids and the Garside normal 
form. For details and proofs we refer to (ECH+92[|Deh02al . 
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In a cancellative monoid M with unit 1, we can define the prefix partial order: 
For x,y G M, we say a; ^ y if there exists c G M such that xc = y. Similarly, we 
define the suffix partial order by saying that x y ii there exists c € M such 
that X — cy. We call s £ M an atom, \i s ^ ah (with a,b G M) implies a = 1 or 
h = 1. We write A for the set of atoms of A/. 

A cancellative monoid M is called a Garside monoid of spherical type, if it is 
a lattice (that is, least common multiples and greatest common divisors exist 
and are unique) with respect to =^ and with respect to if there are no strict 
infinite descending chains with respect to either ^ or ^, and if there exists an 
element A & M, such that D(A) = {s £ Af | s =<; A} = {s G M [ A ;>= s} 
is finite and generates M . In this case, we call A a Garside element and the 
elements of D(A) the simple elements (with respect to A). Moreover, we denote 
the ^-gcd and the ^-Icm oi x,y £ M by x Ay respectively x\/y. It follows that, 
for s G D(A), there exists a unique element ds G D(A) such that sds — A. 

We assume for the rest of this section that M is a Garside monoid of spherical 
type. Since D(A) generates Af , every element x G Af can be written in the 
form X = xi ■ ■ ■ Xm with xi, . . . , Xm G D(A). The representation as a product of 
this form can be made unique by requiring that each simple factor is non-trivial 
and maximal with respect to More precisely, we say that 
in (left) normal form, if Xm 7^ 1 and if Xi — A A (xi • ■ ■ Xm) for i = 1, . . . , m. 
Equivalently, we can require 

(2) Xm 7^ 1 and dxi A Xi+i — 1 for i = 1, . . . , m — 1 . 

If a word is in normal form then all occurrences of A must be at the start, hence 
the normal form of x is of the form A''xiX2 ■ ■ ■ xi where Xi G D(A) \ {1, A}. We 
say that inf(x) = fc is the infimum of x, c\{x) = I is the canonical length of x, 
and sup(x) = fc + Hs the supremum of x. 

As M satisfies the Ore conditions, it embeds into its quotient group Q{M). 
Conjugation by A gives a bijection t: D(A) — D(A) and, as D(A) is finite, this 
implies that some power of A is central. Hence, for every element x of Q{M) 
there exists an integer fc such that A''x lies in M, and so we can extend the 
normal form to the quotient group Q{M). 

Of particular interest to us in [section 3] will be the maps projecting onto the 
i-th non-A factor from the left respectively from the right of the normal form. 
If A''xiX2 ■ ■ ■ xi is in normal form, we define 



Classical Garside structure for the braid group 

The monoid defined by the presentation ([1]) is a Garside monoid of spherical 
type whose quotient group is the braid group i3„ on n strands. It is referred to 
as the classical Garside monoid for i?„. The atoms of are the generators 
ai, . . . , (T„_i, and the Garside element A of i3+ is the so-called half-twist, the 
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positive braid in which any two strands cross exactly once. The simple braids are 
exactly those braids in which any two strands cross at most once. In particular, a 
simple braid is characterised by the permutation which it induces on the strands, 
whence the set D(A) is in bijection to the symmetric group Sn- 

Given x e we define the starting set of x as S{x) = {a E A \ a ^ x} and 
the finishing set of x as F(x) = {a G A \ x )^ a}. 

We remark that the conditions characterising normal forms can be expressed 
in terms of strating and finishing sets, as an element x G is simple if and 
only if it is square-free, that is, if and only if it cannot be written as x = ua^v 
with a,u,v e 5+ and a ^ 1: For any x G D(A), we have F{x) n S{dx) = 
and F{x) U S{dx) = A jCha95[ Lemma 4.2]. For u,v e D(A), one therefore has 
duAv = lii and only if S{v) C A\ S{du) = F{u). 

3 Normal form 

In this section we will investigate the normal form of random elements. For 
our experiments, we constructed and analysed samples of 9999 elements of 
for each combination of number of strands n € {5, 10, 15, 20, 25, 30} and word 
length k e {4,8,12,16,24,32,48,64,96,128,192,256,512,1024,2048} for both 
uniformly random words and uniformly random braids. For uniformly random 
words we also analysed samples with a word length of 4096. The samples of 
uniformly random braids were constructed using an implementation of the al- 
gorithm described in |GGM13] by the first author; the rest of the computations 
were done using a development version of Magma |BCP97) V2.19. 

Using these samples we will investigate the distribution of simple factors along 
the normal form of the elements, that is, we will look at how the induced prob- 
ability measures Ai,(Wordfc), Ai,(URBfc), /9i,(Wordfc) and pi*(URBfc) on the set 
of simple elements vary with i. This will lead us to investigate how the normal 
form changes when an element is multiplied by an atom. 

3.1 Stable region 

The fact that there are a large number of simple elements makes it impractical 
to look directly at the distribution at each position of the normal form. So, we 
will use several invariants instead, namely the word length and the starting and 
finishing sets, to indirectly probe these distributions. 

Word length 

[Figure l] shows how the mean factor length varies along the word. We observe 
that, provided the word is long enough, the word can be divided into three 
regions: An initial region where the word length of the factors is rapidly de- 
creasing; a stable region where the word length is constant; and a terminal 
region where it drops to zero. Moreover, the shape and size of this initial region 
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Figure 1: Mean factor length 
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is independent of the word length. The same structure occurs for the samples 
not shown here. 

The variation in the canonical length within each samples has "smeared out" the 
terminal region, causing it to grow in size as the word length increases. If we 
were to draw right justified plots, that is, if the cc-axis was the distance from the 
end, then you would see that, like the initial regions, the terminal regions have 
a constant size and shape for sufhciently long words. 




I ^ ^ ^ ^ ^ 1 
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Figure 2: Mean factor length inside stable region. 
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Figure 3: Distribution of factor lengths in the stable region, ti = 30 

[Figure 2] shows how the mean factor length inside the stable region depends on 
n for both uniformly random words and uniformly random braids, and [Figure 3| 
shows the distributions of factor lengths in the stable region for n = 30 for both 
uniformly random words and uniformly random braids. 

For uniformly random words, the observed mean factor lengths are consistent 
with a linear function in n (the best fit of a model of the form n'^ is obtained 
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for c w 0.9952), whereas for uniformly random braids the mean factor length 
grows much more slowly; the best fit of a model of the form n'^ is obtained for 
c w 0.44941. 

The data indicates that normal forms of uniformly random words are much 
more "densely packed" than those of uniformly random braids, and that this 
difference becomes more pronounced with increasing n. This is consistent with 
the fact that the distribution of braids obtained by choosing uniformly random 
words is biassed towards multiples of the Garside elements of standard parabolic 
subgroups, that is, the Icms of subsets of A |GGM13) . Such Icms have the 
maximal number representing words, as they can be rewritten using all braid 
relations between the generators involved. On the other hand, these Icms also 
have the maximal possible word length of all simple elements in the standard 
parabolic subgroups, that is, they yield the densest possible packing of the 
involved generators into simple factors. 

Starting and finishing sets 

[Figure 4| shows, for a given generator, the relative frequency with which that 
generator lies in the starting set for each canonical factor. To avoid the problem 
with the variation in canonical length smearing out the end of the words we have 
drawn a left justified plot for the beginning and a right justified plot for the end 
of the word. The plots not shown here for different values of n, different word 
lengths, different generators and for the finishing sets all have a similar shape. 
As we saw for the mean factor length, there is an initial region, a stable region 
and a terminal region. In [Figure 4| for uniformly random braids and generator 
15 there is a local minimum around the 10th factor. Nevertheless, the size and 
shape of the initial, and terminal, regions remains fixed once the word length is 
sufficiently long. Furthermore, the sizes of these regions are consistent with the 
sizes observed for the mean factor length. 

One clear difference between uniformly random words and uniformly random 
braids is in the structure of the starting and finishing sets, see [Figure 5] For 
uniformly random words the distribution is mostly fiat only rising at the edges 
of the braids, whereas for uniformly random braids we have a continuously 
changing distribution which rises from almost zero at the edges of the braid and 
peaks in the middle. 

Combining mean word length and starting / finishing set frequencies 

Given a sample of random braids, we consider the mean factor length and the 
relative frequency of each generator being in the starting set respectively in the 
finishing set as functions / of the position p in the non-A factors of the normal 
form. 

For each of these functions /, we identify the interval [pi,P2] that minimises the 
ratio T^T^^ffjj^, where, for 5 C R, we define \S\ = max(5)-min(5). (Intuitively, 
this procedure locates the "most horizontal part" of the graph of /.) We then fit 
a linear model / to /|[pi,p,] and accept the interval [pi,p2] as stable region for 



7 





25 50 50 25 

Distance from start Distance from end 

(a) Uniformly random words 



0.3 



0.2 



0.1 










Generator 

1 

8 

15 






\ 



0.3 



0.2 



- 0.1 



25 50 50 25 

Distance from start Distance from end 

(b) Uniformly random braids 







Figure 4: Relative frequency of a generator being in the starting set. n = 30, 
word length = 2048. 
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Figure 5: Relative frequency of a generator being in the starting and finishing 
set in the stable region, n = 30 



/ if |/([pi,P2])| < 0.1 • |/([pi,P2])|; otherwise we consider the stable region for / 
as empty. (Intuitively, this procedure ensures that the trend in the restriction 
of / to the interval [pi , p2\ is small compared to the statistical fluctuations of / 
on the interval.) 

The stable region of the sample is taken to be the intersection of the stable 
regions for all the functions. 



[Figure 6] shows the start of the stable region, determined as described above, 
as a function of n and the word length for both uniformly random words and 
uniformly random braids. The data shows that, for fixed n and a given method 
of generating random braids, stable regions exist for sufficiently long random 
braids, and that their starting position does not depend on the word length of 
the random braids. 



Our observations thus lead us to making the following conjecture. 
Conjecture 3.1 (Stable region). Consider the braid monoid for any fixed 
n G N. For Hk = Wordfe, respectively URBfe, and for each i the sequences of 
probability measures Ai^(/ife) and Pi^,{^ik) on the set of simple elements converge 
as fc — > cx). Moreover, there exists a probability measure S on the set of simple 
elements and constants C and D such that 

yi > C Xi^ifik) ^ S as k ^ oo 
\fi > D Pi^i^ik) — i- S as k ^ oo 
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Figure 6: Start of stable region. 
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3.2 Bounded penetration distance 

We can view a uniformly random word as the result of a random process adding 
one letter, chosen at random, at a time. The stable region conjecture suggests 
that the change to the normal form when multiplying by an atom is unlikely to 
penetrate into the stable region. This leads us to investigate how the normal 
form of a random braid changes upon multiplication by an atom. 

The normal form of a word wa, where w is in normal form and a is an atom, 
can be calculated from the normal form of w by working through the word from 
the right to the left, repeatedly applying the rewriting rule xy — >■ {xm){m^^y) 
where m = dxAy. If at any point we have m = 1 then we can stop. If we have 
xm = A then all the following rewrites will be of the form a; A At{x), that 
is, consist of an application of the Garside automorphism; we will consider this 
a trivial change. 

Definition 3.2. For two braids x and y the penetration distance pd(a;, y) for 
the product xy is the number of simple factors at the end of the normal form 
of x which undergo a non-trivial change in the normal form of the product: 

pd(a;, y) = c\{x) - max {i e {0, . . . , cl(x)} : xA" '"^ A A* = xyA' '"^ f"^' A A* } 

Using the same samples of uniformly random words and uniformly random 
braids as before, we took each braid and calculated the penetration distances 
for its product with each generator. The mean penetration distance for each 
word length is shown in [Figure 7[ There are clear patterns here: the mean pen- 
etration distance converges as the word length increases; the value of the mean 
penetration distance increases with n; and it is significantly larger for uniformly 
random braids than it is for uniformly random words. 

Conjecture 3.3 (Bounded expected penetration distance). Consider the braid 
monoid for any fixed n E N. Let be the uniform probability measure on 
the set of atoms. For /ij, = Wordfc, respectively URB^, there exists C such that 
for all k, 

E^.XM^bd] < C 

Corollary 3.4. There exists an algorithm to compute the normal form of a 
braid that has linear expected run time. 

Proof. Consider [Algorithm l| for computing the left normal form of a word. The 
first loop is similar to the usual algorithm jECH+921 Ch. 9] except that: we 
increment a counter / each time a A is created, this will be the infimum of the 
normal form; we add a power of A before each simple element, these are stored 
modulo c the central power of A; and the inner loop stops as soon as a new A 
is created. We then add an additional pass, working backwards along the word, 
to push all the As to the front. 

This algorithm has an invariant: the equality x — A^"^; '•A'ia;iA'2a;2---A''=a;fc 
in Bn, where the summands in the first exponent of A are regarded as integers, 
remains true after each line has been completed. 

As the Ij are elements of ^ the operation of "pushing" the As over a simple on 
lines [9l fTSl and [35l has bounded run time. 
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Figure 7: Mean penetration distance. 
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The inner loop in line [TT] pushes the change to Xi through the prefix of the word 
up to that point, so the contents of the loop will be executed pd(a;ia;2 • • • Xi-i,Xi) 
times. Hence, by |Conjecture 3.3| this loop is expected to take a constant amount 
of time. The two outer loops have at most k iterations with each iteration taking 
a constant expected run time, hence the whole algorithm has linear expected 
run time. □ 



[Figure 8| shows the mean penetration distance of each generator for n = 30 
and a word length of 2048. A similar shape can be seen for the other values 
of n. We see that not only is the mean penetration distance longer for uniformly 
random braids, but also the ratio of longest to shortest is significantly larger: 
For uniformly random words the ratio is less than 2, but for uniformly random 
braids it is greater than 10. 
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Figure 8: Mean penetration distance for each generator, n = 30, word length 
= 2048. 



[Figure 9| shows the distribution of penetrations distances observed in our sample 
for n = 30 and a word length of 2048. 
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Algorithm 1 Calculate the normal form of a word x 
Input: x\X2 - ■ • Xk ^ X where each Xi €: A 
Output: 1^^X1X2 • ■ • Xk 
Zi, Z2, • • • , ^fe e ^ 

/ ^ e z 

i 1 

while i < k do 

if dxi A Xi+i ^ 1 then 

m dxi A ajj+i; aij Xim; Xi+i m~^a;8+i 

if j > 1 then 

a;j_i r'^ (a^^j-i); ^j-i Ij-i + lj\ Ij <r- 
end if 

while Xj 7^ A and j > 1 and dxj^i A .Xj 7^ 1 do 
m ^ dxj-i A Xj-; Xj-i Xj-iw; •(— m^^Xj 

J ^ J - 1 
if j > 1 then 

Xj-i -e- r'-* (xj-i); Ij-i + Z^; Ij <r- 

end if 
end while 
if Xj = A then 

Ij+i ^ Ij + Ij+i + 1; 7 7 + 1 

Delete Xj and moving following terms forward, decreasing by 1. 
« ^ « - 1 
end if 

if Xi+i = 1 then 

Delete Xj+i moving following terms forward, decreasing A; by 1. 
if i = then 

i ^ 1 
end if 
else 

I ^ i + 1 
end if 
else 

i ^ i + 1 
end if 
end while 
for j = fc to 2 do 

Xj-i -S- {Xj-l); Ij-i <r- Ij-i + Ij-, Ij <- 

end for 
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3.3 Garside groups 



Clearly the stable region conjecture and the bounded expected penetration dis- 
tance conjecture make sense in any Garside group and for different sequences of 
probability measures. We will give an example of a Garside group Gi where the 
penetration distance is bounded, in other words there exists a constant C such 
that for any element x and any atom a we have pd(a;, a) < C. This stronger 
condition implies that both the bounded penetration distance conjecture and 
the stable region conjecture hold for Gi. We will then go on to give a method 
to establish whether a variant of the bounded penetration distance conjecture 
holds within a given Garside group. 



A small Garside group 



Let Gi be the Garside group given by the following presentation. 

Gi = {A,B \ ABA = BB) 

This group is isomorphic (as a group) to the braid group on three strands, but 
the Garside structure is distinct [Deh02b| . The Garside element is BBB and there 
are eight simple elements. Figure 10| shows the structure of the left ordering on 
the simple elements. [Figure IT shows the matrix where the xy entry is 1 if 
dx A y = 1 and otherwise. From this matrix one can easily read off which 
pairs of simple element are in left normal form. 
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Figure 11: The entry a, 6 is 1 if 9a A 6 = 1 
and otherwise. 



Figure 10: Hasse diagram 
for the left ordering on the 
lattice of simple elements 

Proposition 3.5. For any element x G Gf in the positive monoid and any 
atom a £ {A,B} we have 

pd{x, a) < 3 

Proof. First consider the atom A and look at how the change upon multiplic- 
ation by A can penetrate through possible suffixes of the left normal form of 

X. 
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If we look at sA for each simple element suffix s we see that there are three 
suffixes with canonical length one that have non-zero penetration distance. 



s 


A 


B 


AB 


BA 


BB 


BAB 


Normal form of sA 


AA 


BA 


BB 


BAA 


BBA 


A 


pd{s,A) 





1 


1 








1 



So for a suffix s with canonical length two to have a penetration distance greater 
than one the last factor cannot be A, BA, or BB. We can also rule out BAB 
as in this case a A is created, which cannot affect any earlier factors. [Figure iT] 
indicates which pairs of simple elements can be adjacent in a left normal form; 
using this information, we can easily produce a list of possible suffixes. 
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So for suffixes with canonical length three there are only two possibilities for 
the last two factors where the penetration distance could be greater than two: 
AAB and BAAB. This gives the following list of possible suffixes. 
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Hence pd(a;, A) cannot be greater than two. 

We will now follow the same procedure for the atom B. 

If we look at each simple element we see that there are four possible suffixes 
with canonical length one that have non-zero penetration distance. 
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So for a suffix with canonical length two to have a penetration distance greater 
than one the last factor cannot be AB, or BAB. We can also rule out BB as a 
A is produced. This gives the following list of possible suffixes. 



s 


AA 


BAA 


BBA 


ABB 


BABB 


ABBA 


BABBA 


N.F. of sB 


AAB 


BAAB 


BBAB 


AA 


ABA 


ABBAB 


BAB BAB 


pd(s,A) 


1 


1 


1 


2 


2 


1 


1 



All words with a penetration distance of two create a A so any preceding factors 
would not be affected. Hence pd(a;, B) cannot be greater than two. □ 

Remark Since the Garside element A is central, any changes to the normal 
form when multiplying by an atom are limited to a fixed number of simple factors 
at the end of the word. This means that the normal form can be computed in 
a single pass with a finite state transducer (an automaton with output) which 
is augmented with an integer counter to count each occurrence of A. 
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Penetration sequences 



In the proof of |Proposition 3.5| we considered sequences of simple elements in 
normal form and how the change penetrated through them when we multiplied 
by an atom. We can formalise this idea as follows. For the rest of this section, 
let be a Garisde monoid with Garside element A. 

Definition 3.6. A word (si, toi)(s2, ■ ■ • (sfc, mfc) e (D(A) x D(A))* is a 
penetration sequence if, for all i, 

(3) s^^l,A ^1,A 

(4) i > 1 =^ s,m^ ^ A 

(5) 5sj A Si+i 1 

(6) TOi = dsi A Si+imi+i 



For a penetration sequence, condition dS]) ensures that siS2 • ■ - Sfc is in normal 
form. If we consider how this normal form changes once we multiple my 
then condition © means that vrii is the simple factor that moves out of Si+i 
into Si. Condition also means that sirfii ^ A. 

We are only interested in the canonical factors of the word and not the initial 
power of A. Also, we are only interested in the region where there is non-trivial 
movement between the factors. So we impose condition ([3]). 

We only allow a A to be created at the very start of the sequence, as otherwise 
the initial terms just correspond to conjugating by A which we consider a trivial 
change. Hence condition Q. 

Let PSeQj, denote the set of all penetration sequences of length k. 

For an element x G and a simple element s e D(A) we will say that a 
penetration sequence (si, toi)(s2, TO2) • ■ • (sfc,mfe) is a ■penetration sequence for 
xs if S1S2 ■ • ■ Sfc is a sufhx of the normal form of x, and = dsk A s. The penet- 
ration distance for the product xs equals the length of the longest penetration 
sequence for xs. 

Let 

<^(fe) ■= {xeG+ : cl{x) = k, ini{x) = O} 
Lemma 3.7. There exist constants a, (3 , p, q > Q such that 

iPSEQfcl e 0{kPa^) 

Proof. Since a word in (D(A) x D(A))* is a penetration sequence if each con- 
secutive pair of letters satisfy conditions (121)-®, we have that IJ^ PSeq„ is 
a regular language. Similarly, IJj, is also a regular language. Hence the 
growth functions of these languages are rational functions. 

By |FS09[ Theorem IV. 9], if /fcz*^ is a rational function with poles pi,p2, ■ • ■ , 
Pm then their exist polynomials Pi, P2, • ■ • , P-m such that fk = Pi{^)PT^ ^ foi' 
k larger than some fixed K. Hence the results hold with a and j3 equal to the 
reciprocal of the radius of convergence of the corresponding growth function. □ 



|G+ |eO(fc^/3'=) 
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Theorem 3.8. Let Vk he the uniform probability measure on G^-^. // the ex- 
ponential growth rate of |G^j| is greater than that of |PSeQj,|, in other words 
(3 > a, and IG^-jj € il[k'^j3'') then the expected value Ei,^x/i^[pd] of the penet- 
ration distance with respect to is bounded. 

Proof. If |PSEQi,| is eventually then the penetration distance is bounded. So 
we may now assume otherwise. 

If a < 1 then for every e > we would have that eventually |PSeQj,| < e. So 
as IPSeq^I takes integer values we have that |PSeQj,| is eventually 0, which is 
a contradiction. So we may now assume that a > 1. 

Let 

Xi,k := {(a;, a) G G^^ x A : pd(a;,a) = i} 

Now 

^ ^ 1^" k I 

Given (a;, a) G Xi^k, let {si,mi){s2,rn2) ■ ■ ■ {si,mi) be a maximal penetration 
sequence for xa. From the definition of Xi^k we have that / = i and, as a is 
an atom, = a. The simple factors si, S2, . . . , Si will be the final i factors of 
the normal form of a;, so a; = xiSiS2 • • • Si for some xi G G^i^_^y This gives an 
injective map Xi_k ^ ^^k-i) ^ PSeq,, therefore 

(7) < |G+„,)|.|PSeq,| 

Bv lLemma 3.71 and the requirement that \G^}~) \ G fl{k'^/3''), there exists a G and 
K such that for all k > K 

iPSEQfel < GFa'= ^fc'^/S'' < |G+)| < Gfc«/3'= 

So ior k — i > K we have 

|G(fe_,)l_[PSEQ,| ^ ^3 (fc - 1)9/3'=-* iPa' 

For k — i < K and k > 2K we have i > K and 



Gffc_,)|-|PSEQ,| i?|PSEQ,| ^^s^J-a^ ^^2n,.p-„/« 



(9) ^ 7 ' , < : 7" < G"i?— 9r < C'DkP-'^ , ^ 

|G+)| - |G+)| - W - 



where D is the largest value of |Gj^ | for j G {1, 2, . . . , K}. 

By making use of (O and splitting the sum into two parts we have that for 
k> 



E..x..[pd]< E ^ I I.I + E 



'=-^'-1 |G +_,)|-|PSeq.| J;^ . |G+_,)|-|PSeq. 

\Gt,,\-\A\ 

« ' V 



18 



By ® 




Which converges as fc ^> oo if a < /3. 
By ®, if fc > 2K then 

Which also converges as fc ^> oo if a < /3. 

Hence, if a < /? then E^,^ x/^^i [pd] is eventually bounded by a convergent sequence 
and so is itself bounded. □ 

Calculating the exponential growth rates 

Let M be the matrix with entries indexed by pairs {s,m) e D(A) x D(A) such 
that s,m ^ {1, A} and sm ^ A whose entries are as follows. 

_ J 1 if S2TO2 7^ A, dsi A S2 ^ 1, mi = dsi A S2m2 
- |q otherwise 

We now have 

IPSeQoI = 1 iPSEQfcl = wM''-\ 

where v is the column vector with all entries 1 and w is the row vector with all 
entries 1. So we can compute the initial terms of the sequence and the minimal 
polynomial of M gives us a recurrence relation. This allows us to compute the 
generating function. The exponential growth rate a is then the reciprocal of the 
radius of convergence of the generating function. 

A similar matrix allow us to compute the generating function and exponential 
growth rate /? for the sequence IG^^j. To verify that \G^^.-^\ G VL{k'i(3^) it is 
sufficient to check that there is a unique pole of the generating function whose 
absolute value is equal to the radius of convergence. 

The results of applying this procedure to several different Garside groups are 
listed in ITable Tl Here Gi — {A, B \ ABA — A^) is the group described above; 
G2 = {x,y,z I xzxy = yzx'^,yzx'^z = zxyzx, zxyzx = xzxyz) is a Garside 
group described in |Pic03j : Bn is the classical Garside structure and BKLn is 
the Birman-Ko-Lee Garside structure for the braid group on n strands. 

Note that in all these examples a < /3 and the generating function for | G^^ | has 
a unique pole with minimal absolute value, hence all of these groups satisfy the 
conditions of [Theorem 3.81 We can also see that, in addition to Gi, the Garside 
groups Bs and BKL3 also have a bounded penetration distance. 

[Figure l"2| shows a plot of a against f3 for these groups. A line is drawn through 
the BKL4 and BKLq points. Although the number of points is small, the fact 
that all of the braid groups lie close to this line is very suggestive. 
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Er=oipsEQ,iz'= 



a 



Gi 

G2 

Bs 

B4 

Bs 
BKL3 
BKL4 
BKL5 
BKLe 



80^ + 130 + 1 



-ST^-" - 102^" 



Z--T 
6z + l 

Az^ - 32z^ + 145/ + 178g ^ 
3 



■780 



lOz'^TlS?"- 7z + 1 



- 333g^ + 97-z + 1 



280z^^ + 



■ + 1592z + 1 



: + l 



3z + l 



Sz' 



52^** + 68^-^ - .56z^ + 23z + 1 



Az" 



- 8?T8? 



-24576000^^ + 



- 5z + 1 
+ 1540 + 1 



16384000^ 360 + 1 

-6640000885716958449490177556480000000000000000000z"^- 
lB16476B48124B908367B1427174400000000000000000000zl''7_ 



-9662-1 



•+209Z-1 





1 



3.532 . 
12.82 . 


3.130. 
8.822 . 
25.31 . 



■ Az ■ 



-llz' 



'Iz - 1 
f 130^ - 



90 



120'' 



160^ + 130 

-20-1 
20 - 1 



60'" - 30-^ 



140 



Qz-' 



150" +80-1 



-1440" + 480^ + 5940'' - 3070^ - 9O0 



1440° 



■ 4800"* + 4980^ 



1990^ + 280 -1 



IO0 



20^ + 40^ 
—3 



•50-= 



+ 40 + 1 
^ 



200'' + 190" -80+1 



-40z^-188z'^+444z'^-482z^ + 122z^-|-104z^-73z'' + 16z+l 
560z» - 1968^7 +3364z<i -3402zS +2132Z* -S36z'^+197z^-2iz+l 



22400000^ 



94O8OOOO0 



TT 



+ 6O0 + 1 



700 + 1 



2 

11.72. 
2 

5.449 . 
18.71 . 

2 

4.839 . 
12.83 . 
35.98 . 



Table 1: Generating functions and exponential growth rates 
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Figure 12: Scatter plot ol the exponential growth rates. 
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